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Relativistic effects and tfie role of heavy meson exchange in deuteron photodisintegration are stud- 
ied systematically for photon energies below the pion production threshold. In a (p/M)-expansion, 
all leading order relativistic one-body and 7r-exchange as well as all static heavy meson exchange cur- 
rents consistent with the Bonn OBEPQ model are included. In addition, one- and two-body boost 
effects have been investigated. Sizeable effects from the various two-body contributions beyond tt- 
exchange have been found in almost every observable considered, i.e., differential cross section and 
■ single polarization observables. 
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I. INTRODUCTION 



Electromagnetic deuteron disintegration is one of the basic processes in order to study various aspects of the strong 
interaction of nucleons in nuclei. For example, a large number of experimental and theoretical papers have clarified 
the role of pion degrees of freedom as manifest in meson exchange currents (MEC) and the importance of relativistic 
' contributions However, the role of heavy meson exchange (with exception of the p meson), which give an important 
contribution to the NN interaction remained largely unclear. One of the reasons lies in the Siegert theorem which 
provides via the Siegert operators in conjunction with the Siegert hypothesis a convenient calculational tool to include 
implicitly the major part of MEC in the electric multipoles 

But admittedly such a procedure overshadows the underlying physics, and only "patches" an inconsistent calcu- 
\ lation. A calculation that uses a set of consistent electromagnetic (e.m.) operators with respect to the hadronic 
interaction model should always be preferred. Such a consistent treatment has been given some time ago allowing 
only pions to interact with nucleons In particular, the importance of consistency of the leading order relativistic 
'"^ I contributions has been pointed out in this work. However, the extension to a realistic interaction model was miss- 
^ ing. Only recently, Levchuk has presented a nonrelativistic calculation of deuteron photodisintegration for the Bonn 
(~| , OBEPR model where all heavy meson exchange currents were included explicitly so that Siegert operators were not 
ILJ ' needed Unfortunately the results are presented in such a way, that the specific contributions from heavy meson 
exchange are not evident. Furthermore relativistic contributions have been neglected, which however become increas- 
ingly important at higher energies. As a sideremark, we would like to mention another earlier calculation based on 
\ the Bonn potential models with a consistent pion exchange current but contributions of heavier mesons were included 
via the Siegert operators only j5||^. 

More recently, we have extended the work of in order to investigate this question for deuteron electrodisintegration 
taking as interaction model the Bonn OBEPQ versions |^,^. As general result we found that the p meson gives 
the most important heavy meson contribution whereas the influence of 77, w, cr, 5, and ^npjio is much smaller, in 
some observables completely negligible, in particular, near the quasifree kinematics. However, it is a priori not clear 
that this conclusion will be valid also for photodisintegration because of the fixed energy-momentum transfer relation. 
Therefore, we want to provide with the present work the missing study of the influence of heavy meson exchange 
on deuteron photodisintegration with inclusion of competing relativistic effects in the one-body and pion exchange 
sector. 

The calculational framework is the same as in |7[| and will be very briefly reviewed in Sect. O. The results are 



presented and discussed in Section III restricting ourselves to the differential cross section dajd£l, and all single 
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polarization observables, i.e., the photon asymmetry S , the target asymmetries T/m, £ {11,20,21,22}, and the 



proton and neutron polarization components Py{p) and Py{n). Section IV gives a short summary and outlook 



II. THEORETICAL FRAMEWORK 

The calculation of the photodisintegration process is based on the equations-of-motion method for the derivation 
of the hadronic interaction model and the corresponding electromagnetic current operators. It is described in detail 
in 1^ where for the first time a consistent treatment including leading order relativistic contributions had been 
presented for a pure one-pion-exchange model. As mentioned above, we have extended this work to the realistic Bonn 
OBEPQ versions and first applied it to electrodisintegration |^ where further details can be found. In particular, all 
explicit expressions for the electromagnetic operators, which can be derived in the equation-of-motion method Q or 
in the unitarily equivalent 5'- matrix approach |l0[| , are listed in the Appendix of [Q, including in addition the jirp- 
and 77rtt'-currents and the currents involving A-isobars. These dissociation and isobar currents introduce additional 
uncertainties and model dependence and will not be considered in the present study. 

When calculating electromagnetic properties of hadronic systems one must, of course, use a set of electromagnetic 
operators that is consistent with the underlying hadronic interaction model, as is demanded by the requirement of 
gauge invariance. In order to construct such a consistent set of operators one should use the same theoretical basis as 
for the hadronic interaction, i.e., for a one-boson-exchange potential (OBEP) one should calculate all corresponding 
meson exchange current (MEC) operators consistently. However, it is not sufficient to simply take the same meson 
coupling constants and cutoffs as used in the potential model. Especially for the pion, care must be taken, because 
there are several sources of unitary freedom in constructing the pionic operators, leading to unitary parameters in the 
corresponding current expressions |]3|,|l0|. These should be chosen consistently with the OBE potential and for this 
reason we will briefiy discuss them now. 

First of all, in view of the unitary equivalence of pseudoscalar (ps) and pseudovector (pv) pion nucleon couplings, 
one introduces a mixing parameter fi, that allows arbitrary mixing of the two coupling types, where /i = corresponds 
to pure ps and /j, = 1 to pure pv coupling. Secondly, it is well known, that this equivalence breaks down in the presence 
of the electromagnetic interaction. Then one has the choice between a ps-pv equivalent theory or one which is not, 
i.e. one can retain or leave out the equivalence breaking term. Clearly, chiral symmetry prefers the pv coupling. Thus 
it is customary to multiply the equivalence breaking term with a new parameter u that switches this term off {v = 0) 
or on {v = 1) and introduces as additional parameter 

7 = /i + i^, (1) 

so that 7=1 corresponds to a chiral invariant interaction theory whereas 7 = violates chiral symmetry. 

Another freedom arises in the p/M expansion for the nonrelativistic reduction, the so-called Barnhill freedom 
described by the Barnhill parameter c [ pT| , which can be incorporated in the definition of a parameter 

/~i = /i + c(l-Ai). (2) 

Note that only for ps coupling one has a dependence on the Barnhill parameter. 

Finally, another unitary parameter stems from retardation. Although the OBEPQ versions are static potentials, 
it is not possible to construct a gauge invariant set of electromagnetic operators that are purely static |10(] . This is 
due to the non-local nature of the e.m. operators when one leaves the nonrelativistic limit, e.g., the charge density 
associated with the pion in fiight fulfills the gauge condition with the other retarded operators Q . We have generated 
the retarded potential through a Taylor expansion of the pion propagator keeping the leading term only 

Vret{k)^Vo{k)A{k^)kl (3) 

where Vo(fc) is the static, nonrelativistic potential, /cq the energy transfer at the vertex, and 

A(fc') = — ^ (4) 

the static meson propagator. Certainly, this approximation is valid only below the pion production threshold. Here 
one has again the freedom to express the energy transfer of the pion by the energy transfers of the individual nucleons 
parametrized by a retardation parameter Vret 
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p - -k^^h^^^ + ^-Vret , (!) (2).2 

«,g _ ftp ftp -1- ^ \^Kq + Kq ) 

= ^(^-Qifc-Q2 + ^^^(fc-(Qi-Q2))'), (5) 

where denotes the energy transfer of nucleon "z" . This freedom can be exploited to ehminate the retarded potential 
in the center-of-mass (cm.) frame by the choice I'ret = \- The retarded e.m. operators must then be constructed 
consistent with this choice. 

With respect to the Bonn OBEPQs, one must be aware that these were constructed from the three dimensional 
Blankenbecler-Sugar reduction of the Bethe-Salpeter equation, thus yielding the nonrelativistic form of the kinetic 
energy operator, while still being a fully relativistic potential. The e.m. operators, constructed within the time-ordered 
perturbation theory or within an equation-of-motion approach |^], contain naturally relativistic one-body operators. 
The connection between potential operators, calculated within these two different approaches, can be made by the 
"Coester" transformation [Q. This leads to the conclusion that the operators in are consistent with the Bonn 
potentials for the choice /i = —1, as was found in 

For the exchange of heavy mesons, the operators given in | [Tc|] are consistent with the Bonn potentials due to the 
simplicity of the corresponding Hamiltonians in leading order. For this reason, no additional unitary parameters 
appear. They should show up in higher order terms, which however, can safely be neglected due to the large meson 
masses. 

For the p MEC we had distinguished in between a "Pauli" current, corresponding to the p contribution propor- 
tional to which can be generated from the tt MEC by substituting terms of the form ct • a by (? x a, and the 
"Dirac" current for the remaining operators. This distinction will also be used later in the discussion. 

With respect to the a meson, we would like to mention that the Bonn potentials OBEPQ (A,B,C) need different 
meson parameters for the isospin T — Q and T — \ channels. Introducing the isospin projection operators, the 
potential (and the e.m. operators) can be viewed as a superposition of an isoscalar and isovector scalar meson for the 
two sets of parameters, yielding effectively four scalar mesons 

= -^{l-n ■f2) V''" + ^ (3 + fi -^2) . (6) 

So, strictly speaking, one has to take into account artificial isovector currents introduced through this procedure, 
which obviously tend to cancel each other. A simpler approximation, i.e., comparing the "pure" a meson exchanges 
of the two parameter sets, also leads to the conclusion that this "anomaly" has no visible effect jj]. 

Finally, for a consistent treatment of leading order relativistic contributions one has to include the wave function 
boost which takes into account the fact that initial and final hadronic states move with different velocities, i.e., their 
rest frames to which the intrinsic motion refers are different. A convenient method for treating this frame dependence 
of the intrinsic motion is to introduce a unitary transformation generated by a boost operator x{P) 

=1 )®e-^'^(^) |P-t ), (7) 



where | pint ) describes the intrinsic wave function in the rest frame |12|. Instead of transforming the wave functions 
themselves, the boost effect is incorporated into the operators by 

e'-^ne-'^ Kn + i[xM^ (8) 

where in the commutator only the nonrelativistic part of 17 has to be considered. 

The operator x can be separated into a kinematic, interaction independent part xo and an interaction dependent 
part XV 

X = Xo + Xv- (9) 

For the two nucleon system one has 

( if.P){p.P) , ^ ^ , ((ai-a2)xp)-P 
^° ^ - [ IQM- + j + 8M^ ' 

whereas a nonvanishing, interaction dependent boost operator exists only for pseudoscalar meson exchange (tt, rf} in 
the case of pseudoscalar coupling 1 13 1^ 



X-. = -(n-ft)^(^)^l-«/^e*'i(P).,-P...?+(l„2,. ,11) 

Therefore, for the Bonn potentials, the potential dependent boost appears with (1 — //) = 2. 
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III. RESULTS AND DISCUSSION 



We have calculated the unpolarized differential cross section and all single polarization observables, i.e., photon 
asymmetry E', the target asymmetries Tn, T20, T21, T22, and the final nuclcon polarization Py for proton and neutron. 
Their formal expressions in terms of the basic T-matrix elements 



are given in |15|, where a denotes the fine structure constant, k the asymptotic relative momentum of the outgoing 
nucleons in the cm. frame, q the photon momentum, E^i the energy of the initial deuteron, and jx the nucleon current 
in a spherical basis. 

For the explicit calculation we have chosen four representative photon energies E-^ = 4.5, 40, 100, and 140 MeV 
covering the region between the maximum of the total cross section and pion production threshold. In order to 
distinguish the different influences from pion, rho, and other heavy meson exchanges, we show their effects in separate 
panels for each observable and each energy. In addition, we show an overview and the potential model dependence 
with respect to the versions A, B, and C of the Bonn OBEPQ. Thus each observable is represented by a figure 
consisting of four columns, one for each energy, and five rows for the overview, tt exchange, p exchange, additional 
heavy meson exchange, and potential model dependence. The notation of the curves is the same for all observables. 

In the first row we present an overview of the following effects: the nonrelativistic one-body current (long-dashed 
curve), the relativistic one-body current (dash-dotted curve), to this added the nonrelativistic tt MEC (dotted curve), 
and finally the total result including all heavy meson exchanges (full curve). 

The next row shows the contributions from tt MEC, starting from the relativistic one-body current (dash-dotted 
curve of the first row) to which first the nonrelativistic tt MEC is added (dotted curve) and then the relativistic tt 
MEC - including retardation corrections - (dashed curve). 

The third row displays the effects of the p MEC: here we start from the relativistic tt MEC (dashed curve of the 
previous row) and include first the Pauli p MEC (short-dashed curve) and then the Dirac p MEC (long-dashed-dotted 
curve). 

The effects of the various heavy meson exchanges are presented in the fourth row. To the relativistic one-body 
current plus relativistic tt MEC and full p MEC (long-dash-dotted curve, as in the third row), we add consecutively 
S MEC (dotted curve), lu MEC (short-dashed curve), a MEC (dashed curve), and finally the 77 MEC (full curve). 

The last row shows the potential dependence of the observable with respect to the different versions of the Bonn 
OBEPQ potential, where the full curve represents the version B, the short-dashed version A, and the dotted version 



Now we will discuss in detail the various observables starting with the differential cross section in Fig. |l|. The 
overview shows that in the maximum of the total cross section, at 4.5 MeV this observable is dominated by the 
nonrelativistic one-body current while only the nonrelativistic tt MEC gives a 10 percent enhancement. Obviously 
relativistic effects and heavy mesons are negligible as well as the potential model dependence. 

At the next higher energy (40 MeV) the nonrelativistic tt MEC becomes comparable to the one-body current. All 
other contributions give a small overall reduction, somewhat more pronounced in forward and backward direction. 
However, if one looks at the separate contributions, one notices a subtle destructive interference of different larger 
effects. First, relativistic tt MEC gives a slight reduction in the maximum but leaves the forward and backward 
directions almost unchanged. Next, from p one sees a strong forward and backward reduction from the Pauli current, 
whereas the Dirac contribution mainly leads to a sizeable enhancement in the maximum which, however, is largely 
cancelled by the additional heavy mesons. Finally, one finds a small model dependence of a few percent. 

Considering now the two higher energies (100 and 140 MeV) one readily notices a dramatic increase of relativistic 
effects. First a sizeable reduction appears from the relativistic one-body current showing the well known effect of 
diminishing the differential cross section at forward and backward angles which comes mainly from the dominant 
spin-orbit (SO) current |p^ . The further reduction from the remaining contributions shows again in detail a strong 
destructive interference. In fact, first the relativistic tt MEC surprisingly enhances the cross section in forward 
direction while then the p Pauli current results in a drastic reduction at both extreme angles. The Dirac contribution 
gives again an overall enhancement but of smaller size. This effect of the Dirac p current is somehow surprising, 
because it is roughly of the same size as that of the Pauli current whereas from the size of the coupling constants one 
would have expected a suppression by a factor of about 50. The additional heavy mesons beyond the p show a much 
smaller influence. Their individual contributions are surprisingly big (up to ~ 5%), in particular compared to their 
role in the parametrization of the NN force and their importance in the electrodisintegration of the deuteron. Most 
prominent is the effect of the 6 meson leading to a reduction of the differential cross section. However, looking at the 
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overall result, the heavy meson exchanges tend to cancel each other. With respect to the potential model dependence, 
one sees now a larger variation, in particular also at forward and backward angles which increases with energy. 

The photon asymmetry in Fig. || is very sensitive to two-body effects, as is long known At 4.5 MeV only 
the nonrelativistic one-body current contributes to and no potential model dependence appears. Then at 40 MeV 
the nonrelativistic tt MEC becomes sizeable as well as the Pauli p current. All other effects, relativistic one-body 
and IT MEC, Dirac p and additional heavy meson effects are very small, as is the potential model dependence. At 
higher energies the relativistic one-body current and the relativistic tt MEC becomes important, too. The first leads 
to a sizeable reduction of the photon asymmetry, the latter to a smaller increase. The p MEC increases the photon 
asymmetry, of which the Pauli current is the most dominant part while the Dirac current is comparably small, 
although its size increases with the photon energy. The influence of the various heavy meson exchanges are much 
more pronounced than in the differential cross section, mainly coming from the 5 MEC. But again the various heavy 
mesons tend to interfere destructively. Also the potential dependence is quite large for 100 MeV and 140 MeV, where 
the OBEPQ version A yields the biggest asymmetry, version B intermediate values, and version C the lowest photon 
asymmetry. Thus one might be tempted to single out one potential against others by comparison with experimental 
data. However, one has to be careful in such a comparison l^], because before drawing definite conclusions as to which 
model should be preferred, one has to study in detail the remaining theoretical uncertainties due to the strength of 
the dissociation and isobar currents. Here, the additional independent measurement of the unpolarized cross section 
and the vector target asymmetry Tn could help in fixing the respective strengths of these contributions. 

The vector target asymmetry Tn, shown in Fig. |[ is extremely small at 4.5 MeV and thus will not be discussed 
furthermore. However, at 40 MeV, Tn develops a strong sensitivity to the relativistic one-body current as well as 
to the nonrelativistic and relativistic tt MEC leading to an overall sign change and a strong increase in absolute size 
around 30 degree. The p MEC shows only small effects, slightly more pronounced for the heavier mesons. Also the 
variation with the potential model remains marginal. For the two higher photon energies one readily sees an even 
more dramatic influence of the various currents. Already the relativistic one-body current shows a strong effect, which 
is partly cancelled from nonrelativistic and relativistic tt MEC. Also visible effects come from the p MEC, mainly 
Pauli, and from the heavy meson sector, again mainly from the 5. The potential dependence is as big as for the 
photon asymmetry Y} . 

The next observable, the tensor target asymmetry T20 in Fig. ^ shows sizeable effects from the various currents 
mainly for the regions around 0° and 180°, except at 4.5 MeV where only a very small influence from the nonrelativistic 
TT MEC is seen. For the higher energies, the largest effect stems from the nonrelativistic tt MEC. In addition, at 100 
and 140 MeV also the relativistic one-body current becomes equally important, particularly strong in the backward 
direction of increasing size. This observable reacts only slightly when the relativistic tt MEC is added, and the same 
is true for the p current. The additional heavy mesons appear more interesting. The lo and a meson exchanges lead 
to a visible reduction at forward and backward angles. Finally, T20 is stable against a change of the potential version. 
This is valid for all tensor target asymmetries. 

The next tensor target asymmetry T21 in Fig. || shows the largest sensitivity to the nonrelativistic tt meson exchange 
at the low energy of 4.5 MeV, although the observable is small in absolute size. Also for the higher energies the only 
relevant contribution besides the one-body part comes from the nonrelativistic tt MEC, but the relative importance 
diminishes with increasing energy. Otherwise, T21 is a very stable observable with respect to all remaining contributions 
and to a potential variation. 

The last tensor target asymmetry T22 in Fig. ^ is like T21 also sensitive to two-body effects of which the nonrelativistic 
TT MEC is dominating again. It leads for photon energies of 40, 100, and 140 MeV to a drastic increase of the 
asymmetry. But this observable shows almost no influence from relativistic contributions, neither from the relativistic 
one-body nor from the relativistic tt MEC. Also the p current shows an increasing contribution with increasing energy, 
mainly from the Pauli part. The additional heavy meson exchange becomes significant, too, with the largest part 
from Lu MEC, leading to a slight overall decrease. The variation with the potential model is negligible. 

Finally, we show the outgoing proton and neutron polarizations Py{p) and Py{n) in Figs. ^ and ^, respectively. 
They show very similar behaviour with respect to the different contributions, except for the two higher energies 
where one notes some differences. At 4.5 MeV the polarization is dominated by the nonrelativistic one-body current. 
There is only a tiny contribution from the nonrelativistic tt MEC, and all other effects are completely negligible. 
The relativistic one-body current shows some increasing influence with increasing energy, but the most important 
contribution is here again the nonrelativistic tt MEC whereas the sensitivity to the relativistic tt MEC is very small. 
The p contribution stems mainly from the Pauli p MEC, leading to a reduction at backward angles. The additional 
heavy meson exchanges lead to a small but increasing effect for higher photon energies. The potential dependence 
becomes sizeable at 100 and 140 MeV, in particular for the proton polarization, and also increases with energy. 

Finally, we would like to mention that boost contributions have turned out to be totally negligible. This is true 
for both the one-body boost operators, see also as well as for the contributions of the two-body boost currents. 
It is in contrast to deuteron electrodisintegration where the one-body boosts cannot be neglected fl^ whereas the 
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two-body boost is small also there 



IV. SUMMARY AND OUTLOOK 



Summarizing our results, we may state that relativistic contributions beyond the leading spin orbit current are 
important and should be included in a consistent manner within a realistic hadronic interaction model. In comparison 
with the work of Rcf. where a pure pionic model had been used, we see that the relativistic effects are roughly 
of the same size. For the differential cross section one notices that the relativistic pionic current in the present work 
has a different angular distribution than in whereas the effects in the photon asymmetry are very similar. These 
differences stem from the different potential models and were expected. 

Of the heavy mesons beyond the pion, the p meson is the most important one confirming earlier investigations. We 
also find that in many cases the Pauli current is dominating its contribution justifying to some extent the neglect of the 
Dirac current. However, we have found some observables for which the Dirac current gives contributions comparable 
to the Pauli current. With respect to the other heavy mesons beyond the p meson, we have found that each individual 
meson shows a sizeable influence on the various observables of deuteron photodisintegration. But taken together, 
we find in many of the observables considered strong destructive interference between them {5, uj, ct, and 77) so that 
then their overall effect is quite small. Nonetheless, there are certain polarization observables, like the tensor target 
asymmetry T20, that obviously demand the incorporation of the heavy meson contributions. 

We have not attempted to make a comparison with experimental data because, first of all, our main interest was 
the study of the relative importance of relativistic effects for a realistic NN interaction model compared to the 
current contributions from heavy meson exchange. Secondly, since at low and medium energies the results for cross 
sections and polarization observables agree essentially with those of other potential models, one can expect the same 
kind of agreement with experimental data for these energies, whereas at higher energies (100 and 140 MeV) isobar 
contributions, particularly from the A resonance, have to be included for a meaningful comparison. Furthermore, 
above pion production threshold, the approximative treatment of retardation in the present work breaks down and 
the correct incorporation of pion retardation in the NN- and iVA-sector becomes important These effects 

are definitively much more crucial than the effects of heavy meson exchange. Also relativistic contributions to the 
excitation of the A resonance itself might show some sizeable influence. These questions will be investigated in future 
work. 
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FIG. 1. The differential cross section for various laboratory photon energies. Notation of the curves; (1) overview: nonrel- 
ativistic one-body current (long-dashed); relativistic one-body current (dash-dotted); nonrelativistic tt MEC added (dotted); 
total result (full); (2) tt meson: relativistic one-body current (dash-dotted); nonrelativistic tt MEC added (dotted); relativis- 
tic TT MEC including retardation (dashed); (3) p meson: relativistic one-body plus complete tt MEC (dashed); Pauli MEC 
(short-dashed); Dirac MEC (long-dashed-dottcd); (4) heavy meson: relativistic one-body current plus complete tt and full p 
MEC (long-dash-dotted); S MEC (dotted); uj MEC (short-dashed); a MEC (dashed); r? MEC (full); (5) potential: OBEPQ 
version B (full); version A (short-dashed); version C (dotted). 
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